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Abstract 

A low-energy non-unitary leptonic mixing matrix is a generic effect of a large 
class of theories accounting for neutrino masses. It is shown how the extra CP- 
odd phases of a general non-unitary matrix allow for sizeable CP-asymmetries in 
channels other than those dominant in the standard unitary case. The za, — > v T 
channel turns out to be an excellent tool to further constrain moduli and phases. 
Furthermore, we clarify the relationship between our approach and the so-called 
"non-standard neutrino interactions" schemes: the sensitivities explored here 
apply as well to such constructions. 
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1 Introduction 



Non-zero neutrino masses are evidence for physics beyond the Standard Model (SM). 
The complete theory accounting for them and encompassing the Standard Model should 
be unitary, as mandated by probability conservation. The light fermionic fields in the 
theory may mix with other degrees of freedom. The complete mixing matrices are 
generically unitary. However, the effective 3x3 submatrices describing the mixing of 
the known fields need not to be unitary. 

We will not particularize to any given model of new physics beyond the SM (BSM) 
and will assume that the total theory of the world is indeed unitary. Nevertheless, 
low-energy non-unitarity may result from BSM physics contributing to neutrino prop- 
agation, when the physical measurements are described solely in terms of SM fields. 

Low-energy non-unitarity may generically appear when the new fields of the BSM 
theory are heavier than the electroweak scale v: this is the framework contemplated in 
this work. Consider a theory with scale M ^> v. The impact of those heavy fields at 
low energies is well described in terms of effective operators of dimension 5 and higher, 
invariant under the SM gauge group, made out of the SM fields active at low energies 
and with coefficients weighted by inverse powers of the scale M. 

The only possible dimension 5 operator of this type is the famous Weinberg oper- 
ator pp, 



(where c^ 5 is the coefficient matrix of 0(1/ M) ), which happens to violate B — L, 
an accidental symmetry of the SM. Upon electroweak symmetry breaking, < >= v, 
this term results in Majorana neutrino masses. Such an operator is characteristic of all 
theories with Majorana neutrino masses, such as for instance the Seesaw model [2J. It is 
very suggestive that the lowest-order effect of high-energy BSM physics may be neutrino 
masses, which indeed constitute -together with dark matter- the first experimental 
evidence of particle physics BSM. 

Such tiny neutrino masses may naturally result through the tree-level exchange 
of heavy particles, which may be either fermions or bosons. The exchange of heavy 
SM singlet fermions is the essence of the minimal Seesaw model (Type I Seesaw) and 
its generalizations. SM fermionic triplets may also mediate light neutrino masses [3]. 
Analogously, the exchange of heavy SM triplet scalars has been as well widely explored, 
as in the Type II Seesaw model and its generalizations [I] . Theories with extra spatial 
dimensions may also induce very small neutrino masses through tree-level exchanges 
of the Kaluza-Klein tower of singlet states [5] • Some supersymmetric frameworks have 
similar mechanisms. Furthermore, radiative models of neutrinos masses forbid tree- 
level exchanges while inducing Majorana masses at the loop level. Which ones among 
all these constructions result at low energies in non-unitary leptonic mixing? 

It is expected in all generality that the tree-level exchange of heavy fermions 
(scalars) will (not) induce low-energy non-unitary contributions. In ref. [7] it was ar- 
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gued that those theories inducing corrections to the low-energy lepton propagators are 
putative sources of a non-unitary leptonic mixing matrix. Specifically, whenever the 
kinetic-energy terms for light leptons acquire flavor- dependent normalizations, wave 
function renormalization is required to recover canonically normalized kinetic energies, 
which often leaves as collateral damage a non-unitary leptonic mixing matrix N, which 
replaces the usual unitary PMNS matrix, Upmns- 

Although the non-unitary mixing induced at low energies in the simplest Seesaw 
model is typically expected to be too small for detection in the near future, this is not 
necessarily true for variants of the Seesaw mechanism [6] , or other theories beyond the 
SM. At the same time, we are about to enter an era of high precision in neutrino physics. 
It is pertinent to ask whether such precision can shed further light on departures from 
unitarity, which would probe the new physics behind. 

In ref. [7] the so-called MUV (minimal unitarity violation) was developed and the 
absolute values of the elements of the matrix N were determined, using data from 
neutrino oscillation experiments and weak decays. It turned out that non-unitary 
contributions, if present, are constrained to be smaller than the percent level. As for 
the phases, a non-unitary mixing 3x3 matrix has three phases in excess over those 
in the unitary case. No information on the size of the phases of the mixing matrix is 
available, neither on the standard "unitary" phases nor on the new non-unitary ones, 
as present oscillation data correspond mainly to disappearance experiments. 

It is the purpose of this work to explore the future sensitivity to the CP-odd phases 
of the leptonic mixing matrix, in case it turns out to be non-unitary. All oscillation 
channels will be analyzed. In particular, it will be shown that CP-asymmetries in the 
Vn — > v T channel are an excellent probe of such new physics. Notice that CP-odd effects 
in that channel are negligible in the standard unitary case, in which the golden channel 
for CP-violation is v e — > v^. On the experimental side, our quest has led us to consider 
several future facilities. Because of the interest of tau detection, measurements at a 
Neutrino Factory [8] will be considered in detail and favored over Super-Beams [9] and 
/3-Beams [TO] , even the highest energy ones. 

Some a priori different avenues for new physics explored in the literature are the 
"non-standard neutrino interactions" or "exotic neutrino interactions" . These are usu- 
ally implemented through the addition of effective four-fermion operators to the SM 
Lagrangian [HI [121 [13] . These operators can affect the production and detection pro- 
cesses or modify the matter effects in the propagation. We will clarify the relationship 
between our framework and those proposals. The channels we will explore and the 
sensitivities we will predict will be shown to generically apply to them as well. 

Section [2] on Formalism attends to its name and also contains the main qualitative 
argument about the observability of effects related to non-unitarity. Sect. [3] explores 
the sensitivity to CP-odd effects induced in the — > v T and the v e — > channels. 
The comparison with the results in "non-standard neutrino couplings" -scenarios is 
performed in Sect. HI In Sect. [5] we conclude. Finally, apendix A introduces a formalism 
to derive oscillation probabilities in matter with constant density, in the MUV scheme. 
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2 Formalism 



Let us parameterize the general non-unitary matrix N, which relates flavor and mass 
fields^ 

u a = N ai Vi , (2) 
as the product of an hermitian and a unitary matrix, defined by 

N=(1 + V )U, (3) 

with rj' = rj. The bounds derived in Ref. [7] for the modulus of the elements of NN^ 
from universality tests, rare lepton decays and the invisible width of the Z, also apply 
to the elements of rj, since NN^ = (1 + rf) 2 « 1 + 2rj and it follows that 



(4) 



at the 90% confidence level. The bound on rj^ has been updated with the latest 
experimental bound on r — > /17 [H]. Eq. (j2J) shows that the matrix N is constrained 
to be unitary, within a 10 -2 accuracy or better. The unitary matrix U in Eq. ([3]) can 
thus be identified with the usual unitary mixing matrix U = Upmns, within the same 
accuracy. The flavor eigenstates can then be conveniently expressed as^| 
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It follows that the neutrino oscillation amplitude, neglecting terms quadratic in 77, is 
given simply by 

< u p \v a (L) >= Aif{L) (1 - Vaa - m ) + Y, (v^A^iL) + , (6) 

with 

Alf{L) 4 M \v s a M {L) > (7) 

being the usual oscillation amplitude of the unitary analysis. 

New CP-violation signals arising from the new phases in rj, require to contemplate 
appearance channels, a 7^ (3. The best sensitivities to such phases will be achieved in 
a regime where the first term in Eq. ([6]) is suppressed. This happens at short enough 
baselines, where the standard appearance amplitudes become vanishingly small while 
Aaa(L) — 1. The total amplitude is then well approached by 

< u p \u a (L) >= A S J(L) + 2r)* p + 0(riA), (8) 



1 Throughout the paper, Greek (Latin) indices label the flavor (mass) basis. 

2 As neutrino masses are forbidden in the SM, the handy superscript SM is an abuse of language, 
that we allow ourselves to describe the flavor eigenstates of the standard unitary analyzes. 



3 



where 0(rjA) only includes appearance amplitudes and rj components with flavor 
indices other than a(3. This is an interesting property as it implies that, at short 
enough baselines, each oscillation probability in a given flavor channel, P a p, is most 
sensitive to the corresponding r] a p. The other elements of the rj matrix can be safely 
disregarded in the analyzes below, without implying to assume zero values for them. 
That is, their effect is generically subdominant, a fact that will be numerically checked 
for the main contributions, as explained later on. This also means that the subleading 
corrections from the cross sections and fluxes discussed in Ref. [TJ, which induce also 
0(j] A) corrections, do not need to be taken into account. 

For instance, in a two family scenario and within the above-described approxima- 
tion, the oscillation probability would read: 



P af3 = sm 2 (20)sin^ I —J - 4|r/ a/3 | sm5 af3 sin(2fl) sin I— +4|^|', (9) 

where A = Am 2 /2E and r) a p = \r} a f3\e~ zSal3 ■ The first term in Eq. §§§ is the usual 
oscillation probability when the mixing matrix is unitary. The third term is the zero- 
distance effect stemming from the non-orthogonality of the flavor eigenstates. Finally, 
the second term is the CP-violating interference between the other two. Notice that 
the latter is linearly sensitive to both phases and moduli, a fact which will be at the 
origin of the improvement in the sensitivity to the moduli, for non-trivial values of the 
phases 5 al 3, as compared to previous analyzes in the literature. In this same two-family 
scenario, a CP-asymmetry can be written as 



a cp _ Pa P ~ P &P -MVapl si n <W / 1n \ 
^ P a ? + P^ ~ sin(2fl) sin (^) ' 



where, for illustration, it is implicitly assumed that we work in a regime in which the 
term quadratic in rj in Eq. ([9} is negligible with respect to the first term, that is, with 
respect to the standard contributions. 

In the sections below, we will analyze the new sources of CP violation in the u e — ► 
v T and the — ► v T appearance channels, since present constraints on r\ eil are too 
strong to allow a signal in the v e — > channel (see Eq. (141)). When numerically 
computing a given P a p, the only approximation performed will be to neglect all rj 
elements but that corresponding to the channel under consideration, rj a p. They should 
be indeed subdominant, as illustrated by Eq. (|S]). Furthermore, we have checked this 
approximation as follows. The numerical fits have been performed in two ways. First, 
setting to zero all the elements of r\ except rj a p. Next, allowing the other off-diagonal 
elements of rj to vary (except for rj efl , which is extremely well constrained). The results 
are indistinguishable within the accuracy explored. The effect of the diagonal elements 
is expected to be even smaller since they cannot induce CP asymmetries. 

Finally, the parameterization of iV in Eq. ([3]) is on purpose very similar to that used 
to study non-standard neutrino interactions from four-fermion operators in Refs. [TTl 
Some CP-odd effects have also been considered in these scenarios [T2] • We leave 
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to Section 0] the task of clarifying the similarities and differences with our framework, 
their equivalence in certain regimes and the range of application of our numerical 
bounds to those constructions. 

3 Sensitivity to the new CP-odd phases 

As suggested by Eqs. (ED) and (flOl) . the best sensitivities to CP- violation will be achieved 
at short baselines and high energies, where the standard term is suppressed by sin 2 (4r). 
We will therefore study a Neutrino Factory beam resulting from the decay of 50 GeV 
muons, to be detected at a 130 Km baseline, which matches for example the CERN- 
Frejus distance. For these values, sin( A3 2 lL ) — 1.7- 10~ 2 and sin( A2 2 lL ) ~ 6- 10~ 4 , where 
Ajk = (m 2 — m\)/2E. All terms in the oscillation probability Eq. can then be of 
similar order for the channels z/ M — > v T and v e — > v T) if the corresponding 7] a p values are 
close to their experimental limits in Eq. (J4]). In what follows, we will assume 2 • 10 20 
useful decays per year and five years running with each polarity. 

The appearance of v T s will be contemplated assuming a 5 Kt Opera-like detector. 
The efficiencies and backgrounds for the measurement of v e — > v T transitions at a 
Neutrino Factory have been taken from Ref. [T5|. A similarly detailed analysis for the 

— > v T channel is still lacking in the literature. In fact, — > v T oscillations are 
one of the main backgrounds in the detection of v e — > v T transitions. As the signal 
for the former transition is stronger, looser cuts should be sufficient when analyzing 
it. Namely, since charge identification should not be mandatory, it may be possible to 
study not only the r into /i decay mode, but also the rest of the r decay modes, gaining 
a factor of 5 in sensitivity which, on the other side, will mean larger backgrounds. We 
have thus considered sensitivities and backgrounds a factor 5 larger when analyzing 
the — > v T channel than those used for the v e — > v T channel [T6] . 

In the numerical analysis, the complete oscillation probabilities in matter have 
been used (albeit with the simplifications on rj previously described). It is inter- 
esting, though, to understand qualitatively and in detail the role of matter effects 
and, whenever relevant, the dependence on small neutrino parameters, such as A 21 L 
and sin2# 13 . For this purpose, we will consider below an expansion of the oscillation 
probabilities at higher order than that implied by Eq. ([9]). Taking into account that 
A31L ~ AL ~ 10~ 2 , where A = ^/2Gpn e , with Gp being the Fermi constant and n e the 
electron density in the Earth crust, A 21 L ~ 10^ 3 5 , sin26'i3 < 10~ a5 , \r) a p\ < 10 -2 and 
\r] eiM \ < 3.6 • 10~ 5 , it is consistent to expand to second order in the following parameters: 

sin 2 2^3, A 21 L, (A 31 L) 2 , (AL) 2 , Va/3 , (11) 

with rj efl set to zero, as its contributions are always suppressed by extra small param- 
eters. 
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Figure 1: Left: 3a contours for two input values of |^ MT | and 5^ T represented by the 
stars. Right: the solid line represents the 3a sensitivity to |77 Mr | as a function of 5^ T , the 
dotted line the 3a sensitivity to 5^ T and the dotted-dashed line represents the present 
bound from r — *■ fi'y. 



3.1 The — ► v T channel 



The expression for P^ T expanded to the order just described can be found in |Appendix A 
Matter effects are subleading, as this probability is not suppressed by small standard 
parameters such as sin# 13 or A 12 . The two family approximation in Eq. (Q is thus 
very accurate to understand qualitatively the results and reads 

P^r = sin 2 (2# 23 ) sin 2 f^f^J ~ 2|^ T | sin V sin(2fl 23 ) sin (A 31 L) + A\n^ T \ 2 . (12) 

This equation indicates that the CP-odd interference term is only suppressed linearly in 
\rjfj_rl. This can indeed be observed in the result of the complete numerical computation, 
Fig.[[] which shows the sensitivities to \i]^ T \ and <5 MT obtained. The left panel represents 
two fits to two different input values of \r]^ T \ and 5^ T (depicted by stars). The dashed 
lines correspond to fits done assuming the wrong hierarchy, that is the opposite sign 
for A31 to that with which the number of events were generated. As expected from 
Eq. (TI21) . a change of sign for the mass difference can be traded by a change of sign for 
5 n T . Nevertheless, this does not spoil the potential for the discovery of CP violation, 
since a non-trivial value for is enough to indicate CP violation. Furthermore, 
the sinusoidal dependence implies as well a degeneracy between 5^ T — > 180° — 5 Mr , as 
reflected in the figure. 

The right panel in Fig. [T] depicts the 3a sensitivities to \rj^ T \ (solid line) and 5 Mr 
(dotted line), while the present bound from r — > /ry is also shown (dashed line). 
The poorest sensitivity to [rj^rl, around 10 -3 , is found in the vicinity of 5 MT = and 
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<5 MT = 180°, where the CP-odd interference term vanishes and the bound is placed 
through the subleading \f]^ T \ 2 term. The latter is also present at zero distance and its 
effects were already considered in Ref. [7], obtaining a bound of similar magnitude. 
The sensitivity to \ti^ t \ peaks around \f]^, T \ — 4 ■ 10 -4 for 5^ T ~ ±90°, where sin^ 
is maximum. That is, for non-trivial values of 5^ T not only CP-violation could be 
discovered, but values of |^ T | an order of magnitude smaller could be probed. 



3.2 The u e — > v T channel 

Contrary to P MT , all terms in P eT , standard and new ones, are suppressed by at least 
one of the two small standard parameters sin #13 and A12. The oscillation probability, 
expanded to second order in the parameters in Eq. (jTTi) . reads 



P eT = c 2 23 sin 2 29 i3 ( + s 2 23 sin 2 26 12 1 21 



2 J z * V 2 

A 21 L\ /A31L 



C13 sin 26*12 sin 26*23 sin 26*13 cos 8 ( , } 

- 2|?7 eT |c23 sin 2^i 3 sin(5 + <5 eT )(A 3 iL) 
+ 2|?7 eT |s23Ci3 sin 26 12 sin 5 er sin(A 2 iL) 

+ MVer\\ (13) 



where s^-, c^- stand for cos%, sin%, respectively. 

The first three terms in Eq. (fT3"l) are the usual unitary contributions, suppressed 
quadratically in sin 9 is and A12. The next two are interference terms between the 
unitary oscillation contribution and the non-unitarity parameter r] eT , and are CP-odd 
and suppressed H by \r] er \ and either sin 6*13 or Ai 2 . The last term is the zero distance 
effect only proportional to \r] er \ 2 . Notice that Eq. (jT3l) would also be valid for the 
P eAt oscillation probability replacing s 2 3 — * — c 2 3, c 2 3 — > s 2 3 and r\ er — > r] efl . Although 
in the numerical analysis we have used the full oscillation probability in matter, the 
approximation in which Eq. f|T3l) has been obtained shows no sensitivity to matter 
effects. 

Fig. [2] shows the sensitivities to \r] er \ and 5 er with this setup. The left panel depicts 
a fit to the input value represented by the star: the dependence on 5 eT is seen to be very 
mild and no measurement of this quantity can be performed. This is easily understood 
from Eq. (fT3"j) . for \rj eT \ > 10~ 3 : the last term, proportional to |^ e r| 2 , dominates over 
the CP-violating ones and no information on 5 eT can be extracted. In consequence, 
in the right panel we find 3a sensitivities to \i] eT \ around 10 -3 , but no sensitivity to 
5 eT . The sensitivity to \r] eT \ from the \r) eT \ 2 term is also present at zero distance and 
already studied in Ref. [7j, where a similar bound was obtained. We have checked that 

3 This dependence was previously observed for v e — > v„ transitions in a related context in Ref. |12j . 
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Figure 2: Left: 3cr contour for an input value of \r] eT \ and 5 eT represented by the star; 
the \rj eT \ 2 term dominates and there is no sensitivity to the CP-violating phase. Right: 
3a sensitivity to \r) eT \ as a function of 5 er , the dashed line represents the present bound 
from r — > ej. 



increasing the statistics by a factor 100, values of \r] eT \ an order of magnitude smaller 
would be accessible and the CP-violating terms would start to dominate, providing 
sensitivity to 5 eT . 



4 Non-unitarity vs non-standard interactions 



Non-standard neutrino interactions are usually introduced through the addition of 
effective four-fermion operators to the SM Lagrangian [Til [T2l [13], made out of the 
SM fields and invariant under the SM gauge group. They can affect the production or 
detection processes or modify the matter effects in the propagation, depending on the 
operator or combination of operators considered. 

In order to clarify the impact of each operator, let us first write the effective SM 
Lagrangian after integrating out the W field @. In the flavor basis, it reads 



-C 



cc 



4G F 

7f 



[J2( P a M ^PLla)(d^P L U P ) + {v S a M l,P L l a ){hl,PL4 M ) + ^ ■ ( 14 ) 



aii 



The first term in Eq. (f!4l) contributes to hadronic detection and production processes 
such as v a d — > l~u. The second term describes leptonic production and detection 



4 For simplicity, we do not show the neutral current interactions nor quark sector, since they are 
not relevant for this discussion. 
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processes involving two neutrinos, such as muon decay at the NF, /i + — > e + h , e h , fl , as 
well as matter effects for a = (3 = e. 

Non-standard neutrino interactions are usually introduced through analogous ef- 
fective operators that can be non-diagonal in flavor. For instance, in Ref. [T3J new 
matter effects stemming from the operator ~ (Vf Mr )^ l PL l 'p M ) (flfif) were studied. 
This kind of operators give rise to non-diagonal corrections to the effective potential 
for the evolution equation in matter (in the SM flavor basis). 

In Ref. [12], another kind of four-fermion operators were introduced: 



P = ^E^^^K^M+^.c. 



D = ^GJ(^)(Vl«)U.c., (15) 



which affect the production (O p ) and detection (0 D ) processes at a Neutrino Factory, 
but do not correct matter effects. In this case, the relevant terms in the interaction 
Lagrangian are modified as follows: 



- = ^ E i <W + ~£ I (d^P L u) + h.c. 





+ U«e + ^ (^ M l,P L e) (m,P L ^ M ) + h.C. 



(16) 



Defining now = U a i z/j, with U being the PMNS matrix, the effective production 
and detection states are given by: 

K> = {\ + e^)e P \4 M >=^ + ^ p Uu; l \u l > 

K> = (l + e*Vl^ M >=( 1 + e *V^I^>> (17) 

qP(D) 

where e^f 1 = — , up to normalization factors. 

These expressions are very similar to our parameterization of the effects of a non- 
unitary mixing matrix, Eq. ([5]). In fact, the latter can also be encoded in terms 
of effective four-fermion operators, after integrating out the W and Z bosons. The 
difference is that, in the case of a non-unitary mixing matrix, the coefficients of the 
different effective operators induced by it and contributing to production, detection and 
matter effects are not independent but related. For instance, it follows from Eq. (J3]) 
that ri Ql3 = r]p a . Which would mean e p a g = e^* a , a constraint usually not required 
when introducing non-standard interactions. When such equality holds, the oscillation 
physics induced by non-standard interactions is equivalent in vacuum to that stemming 
from non-unitarity in the MUV scheme. Furthermore, even if no such relations among 
e^g and are assumed, the order of magnitude of the bounds obtained in the previous 
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chapters should apply as well to non-standard interactions, barring very fine tuned 
cancellations. The new CP signals analyzed in this work are also probes of the phases 
of non-standard interactions. 

Finally, in Ref. [7j it was shown that matter effects are also modified in a very 
definite way in the presence of a non-unitarity mixing matrix. These new matter effects 
are not independent from the production/detection effects, contrary to the customary 
assumption in studies of non-standard neutrino interactions. Nevertheless, as we have 
studied a setup such that matter effects - standard and new ones - are negligible, the 
bounds derived here can be applied directly to the case of non-standard interactions. 

5 Conclusions 

It is important to analyze low-energy data without assuming a unitary leptonic mixing 
matrix, as it is a generic window on new physics, even if the effects are expected to be 
extremely tiny in the simplest models of neutrino masses. 

A non-unitary matrix has not only more moduli than a unitary one, but also sup- 
plementary phases which may lead to new signals of CP-violation. In particular, an 
asymmetry between the strength of — > v T oscillations versus that for — > v T has 
been shown to be a beautiful and excellent probe of new physics, when measured at 
short-baselines (~ 100 km.) using a Neutrino Factory beam of energy O(20GeV). 
Non-trivial values of the new phases can lead not only to the discovery of CP- violation 
associated to the new physics, but also allow to probe the absolute value of the moduli 
down to 10~ 4 . This means an improvement of an order of magnitude over previous 
analyzes of future facilities, which were sensitive only to the moduli. 

Our analyzes for future sensitivities to non-unitarity, as well as the new signals of 
CP-violation explored here, have been shown to also apply to the physics of "non- 
standard neutrino interactions" , barring extremely fine-tuned cancellations. 

Appendix A P^ T expansion 

In this Appendix, a convenient formalism is introduced to derive oscillation probabili- 
ties P a p in matter with constant density, in the MUV scheme, and the results are then 
applied to the \x — > r channel. 

Several years ago Kimura, Takamura and Yokomakura derived a nice compact for- 
mula [T7j for the neutrino oscillation probability in matter with constant density. They 
showed that the quantity U*jUpj in matter can be expressed as a linear combination 
of the quantity U*jUpj in vacuum, where U a j and U a j stand for the matrix element 
of the PMNS matrix in vacuum and in matter, respectively. A simple derivation of 
their formula and its generalization are given in Ref. [18J. Here it is shown that the 
framework in Ref. [18] can be applied also to the MUV case. 
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Time evolution of the neutrino mass eigenstate \l/ m = |i/a)j \ v s)) T in matter 

with constant density in the MUV framework is given by 



dt 



(£ + N AN*) ty m (t), 



where £ = diag(£'i, E 2 , E 3 ) is the energy matrix in the mass eigenstate in vacuum, 
A = \/2GFdiag(n e — n n /2, — n n /2, — n n /2), n n is the neutron density, and the non- 
unitary matrix N relates the flavor fields to the mass fields (cf. Eq. (J2j)). Diagonalizing 
the hermitian matrix £ + N T AN* with a unitary matrix W 

£ + N T AN* = W£W~\ (19) 

where £ = diag(i?i, E 2 , E 3 ) is the energy matrix in matter in the MUV scheme, the 
mass eigenstate at distance L can be solved as 

* m (L) = Wexpi-it^W-^UO). (20) 

We define the modified amplitude A(y a -> v p ) = A{y a -> v p )(NN^) l £{NN^ , 

A(v a -> ^) = [j\rW , exp(-i£L)W ,_1 iV T ] a/3 , 

and the modified probability -P(f Q — > = \A(u a — > z^a)] 2 is given by 

P(u a ^up) = |(iV*iV T )^| 2 -4^Re(XfX^)sin 2 (AE, fe L/2) 

+2 £ Im(Xf Xf *) sin(A£ ife L), (21) 

i<fc 

where AE jk = E d - E k and Xf = (N*W) aj {NW*)p j {j = 1,2,3). As in the case of 
the standard neutrino scenario [13 [18], Xj 13 can be expressed in terms of the quantity 
■gap = U*jUi3j in vacuum and Ej. To show it, let us first note the following relations: 

Y^Ej) m Xf = Y,(N*W) aj (E j ) m (NW*) /3j = [N*(£ + N T AN* ) m N T ] a p = Y r 
j j 

form = 0,1,2. (22) 



a(3 
m+1 



Secondly we rewrite Eqs. (122]) as 

3 

£ V jm X^ = Yf for j = 1, 2, 3, (23) 



m=l 



where Vj m = (E m y 1 is the element of the van der Monde matrix V. The simultaneous 
equation (123~1) can be easily solved by inverting V: 



m=l 



(1/AE 21 AE 31 )(E 2 E 3 , 


-(E 2 + E 3 ), 






{l/AE 21 AE 32 ){E 3 E x , 


-(E 3 + E x ), 






{1/AE 3X AE 32 ){E X E 2 , 


-(E x + E 2 ), 


1) / 


V n Q/3 / 



(24) 
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In the present paper we focus on the short distance case, i.e., the case where 
\AEjkL\ <C 1, \AEjkL\ <C 1, \Gpn e ^ n L\ <C 1. It turns out that terms of order (AEj k L) 3 
or higher are negligible at short distance and we expand the sine functions in the 
probability (12"T1) to quadratic order in the argument AEj k L: 

sm 2 (AE jk L/2) ~ (AE jk L/2) 2 , sin(AE jk L) ~ AE jk L 

Plugging Eq. (|24p in the probability (|2ip and using the form of in Eq. (1241) . we 
can show that the first terms take the following simple forms: 

^Re{Xfxf*)(AE jk L/2f = {L/2) 2 [Re{YfYf*)-\Yf\\ 

j<k 

^\m{Xfxf*)AE jk L = -Lhxi{YfYf*) 

j<k 

Hence we obtain the probability at short distance 

P(u a -> up) ~ |(iViVt) Q/3 | 2 - L 2 [Re(Y 1 Q/3 Y 3 a/3 *) - |Kf| 2 ] - 2Llm(YfYf*). (25) 

We are interested in the MUV scheme with small values of rj, where rj was defined 
in Eq. ([3]), so we evaluate Y^ only to first order in 77: 

Yf ~ 5^ + 2^, 

3 3 

j=2 7 j=2 

3 3 

F 3 a Vo = ^(A,!) 2 ^ + ^A il Xf (^ a + ^) + (X / ,) 2 , (26) 

J=2 j=2 

where we have used the energy matrix £ — Eil + N T AN* instead of fTIT?]) for ease of 
calculation, since the shift (—Eil) only changes the phase of the amplitude A{y a — > up). 
Kj was evaluated only in the limit T] a p — > 0, as it always appears together with 3^ , 
which is first order in r) a p. 

It is straightforward to obtain the components Y^ T for j = 1,2,3 from Eqs. (1261) . 
Keeping only terms to second order in sin 2 0i 3 , A 21 L, (A 31 L) 2 , (AL) 2 and first in r) a p, 
and setting r\ eil = 0, we get the probability P^ T = P{u^ — > u T )/(NN^) tltl (NN^ ! ) TT : 

P^ T = sin 2 2^3(A 3 ii:/2) 2 -2|^ T |sinVsin2^3(A3iL)+4|^ r | 2 

- (l/2)c 2 2 sin 2 2MA 3 iL)(A 21 L) " \Vf,r\sm26 23 cos5^(AL)(A 31 L) 
+ (1/4)4 sin 2 2^ 23 (A 21 L) 2 + 2|^ r |c 2 2 sin V sin 20 23 (A 21 L) 

- (l/4)s 13 sin 4# 23 sin 26 12 cos 5(A 31 L)(A 21 L). (27) 

Notice that matter effects in this channel are clearly subdominant. All terms but those 
in the first line of Eq. (|27p are numerically smaller than 10~ 4 for the setup considered 
in this work and have thus negligible practical effects. 

For the e — > r channel the same formalism has been used to obtain the expanded 
probability, presented in Eq. (TIB"]) . 
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